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Abstract 

In this paper we report on our theoretical studies of various aspects of the 
correlated stopping power of two point-like ions (a dicluster) moving in close 
but variable vicinity of each other in some metallic target materials the lat- 
ter being modelled by a degenerate electron gas with appropriate densities. 
Within the linear response theory we have made a comprehensive investiga- 
tion of correlated stopping power, vicinage function and related quantities for 
a diproton cluster in two metallic targets, aluminum and copper, and present 
detailed and comparative results for three approximations to the electron gas 
dielectric function, namely the plasmon-pole approximation without and with 
dispersion as well as with the random phase approximation. The results are 
also compared, wherever applicable, with those for an individual projectile. 
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1 INTRODUCTION 



The stopping of energetic charged particles in a target material is a prob- 
lem of long-standing theoretical and experimental interest. Early pioneering 
theoretical work by Bohr who had a lifelong interest in this problem was fol- 
lowed by Bcthe and others. There is by now an extensive literature on this 
topic. We refer the reader to two recent review articles [1, 2]. The problem in 
its generality is rather complex. However simplified theoretical models have 
been applied with considerable success in explaining experimental data. As- 
suming a weak coupling between the energetic particle and a target material, 
specially a metal which is usually approximated by a degenerate electron gas, 
a detailed theoretical model has emerged through the works of Bohm, Pines, 
Lindhard, Ritchie and other authors [3-5]. 

A comprehensive treatment of the quantities related to inelastic particle- 
solid and particle-plasma interactions, e.g. scattering rates and differential 
and total mean free paths and energy losses, can be formulated in terms of 
the dielectric response function obtained from the electron gas model. The 
results have important applications in astrophysics [6, 7] and radiation and 
solid-state physics [8-11], and more recently, in studies of energy deposition 
by ion beams in plasma fusion targets [12-16]. 

One can think of several situations in which the projectile beam ions 
may be closely spaced so that their stopping is influenced by their mutual 
interactions [14, 15] and thus differs from the stopping of charged particles 
whose dynamics is independent of each other. This can happen, for example, 
in the case of very high density beams, or more realistically, when cluster 
ions are to be used instead of standard ion beams [16]. 

The stopping of uncorrclated or independent charged particles in a de- 
generate electron gas (DEG) has been extensively studied in the literature 
(see, for example, [3-9, 13] and other references therein). These studies have 
been done mostly within the linear response formalism and for the projectile 
velocity V comparable or greater than vp, the electron Fermi velocity. 

The objective of this paper is to make a study of the stopping power 
(SP) of correlated charged particles in a DEG. The simplest and yet phys- 
ically relevant case is the SP of an ion pair (a dicluster). We report on a 
comprehensive investigation, which is mostly numerical, of various aspects 
of a dicluster stopping in a DEG. Previous theoretical works have consid- 
ered this problem within the linear response theory [9-11] and in a simple 
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plasmon-pole approximation. In our study, again in the linear response for- 
malism, we have used both the plasmon-pole approximation (PPA) as well 
as the full Lindhard expression for the random phase approximation (RPA) 
in the DEG dielectric function. As in earlier studies we consider a dipro- 
ton cluster as a projectile and compare our theoretical results with those of 
Basbas and Ritchie [10] who used PPA and with those obtained in RPA. 
Whenever applicable, the results arc also compared with those for an inde- 
pendent projectile (e.g. those due to Yakolev and Kotel'nikov [7]). No RPA 
results for a diproton cluster SP and related aspects have previously been 
reported in the literature. 



2 STOPPING POWER 

Let us consider an external charge with distribution Pext(r, ^) = Qext(r — 
Vt) moving with velocity V in a medium characterized by the longitudinal 
dielectric function e{k. u). Within the linear response theory and in the Born 
approximation the scalar electric potential (p{r, t) due to this external charge 
screened by the medium is given by [1] 

where (5ext(k) is the Fourier transform of function (5ext(r)- 

The stopping power which is the energy loss of the external charge re- 
garded as a projectile, per unit path length in the medium regarded as a 
target material, can be calculated from the force acting on the charge. The 
latter is related to the induced electric field Ei^d in the medium. For a three- 
dimensional medium we have, for the SP, 



S=-J drgext(r-Vt)^Ei,d(r,t) = ^^ j dk |gext(k)|' ^Im-^^i^ 

' (2) 

Eq. (2) is applicable to any external charge distribution. We shall apply 
it to a cluster of two point-like ions having charges Zie and separated 
by a variable distance R. For this dicluster 
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|Qext(k)f = — [Zl + Zi + 2Z,Z2CosikR)] . 



(3) 



Then the stopping power of a dicluster may be written as 



S - {Z^ + Zl) S,^a{V) + 2ZiZ2>5eorr(R, V), 



(4) 



where Si-ad{V) and 5'corr(R', V") stand for individual and correlated SP, respec- 
tively. Prom Eqs. (2) and (3) 



Jq{x) is the Bessel function of the first kind and zero order and ■& is the angle 
between the interionic separation vector R and the velocity vector V. Eqs. 
(5) and (6) can also be obtained from the linearized Vlasov-Poisson equations 
for a two- ion projectile system, as has been done by Avanzo et al. [14]. In 
their study the target material is a dense classical electron gas. 

We note that there are two contributions to the SP of a two-ion cluster. 
The first one, given by the first term in Eq. (4), is the uncorrelated particle 
contribution and represents the energy loss of the individual projectiles due 
to their coupling with the target electron gas. The second contribution, the 
second term in Eq. (4), arises due to a correlated motion of the two ions 
through a resonant interaction with the excitations of the electron gas. Both 
terms are responsible for an irreversible energy transfer from the two-ion 
projectile system to the target electron gas. 

In many experimental situations, clusters arc formed with random orien- 
tations of R. A correlated stopping power appropriate to this situation may 
be obtained by carrying out a spherical average over R of the 5'corr(R, V) in 
Eq. (4). We find 
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(7) 



One may consider an interference or vicinage function, which is a measure 
of the separation of single-particle contribution from its correlated counter- 
part, to the stopping power. This function is defined as [10, 14, 15] 



r(R,y) = 

Eq. (4) can then be put in the form 

S = {Zl + Zl) 5ind(V) 



'5'corr(R-, y) 
^ind(^) ' 



i + .^r(R,y) 



(8) 



(9) 



Zl + Zl 

r(R, V) describes the strength of correlation effects with respect to an un- 
correlated situation. The vicinage function becomes equal to unity as 7? ^ 
when the two ions coalesce into a single entity, and goes to zero as — > oo 
when the two ions are totally uncorrelated. 



3 THEORETICAL CALCULATIONS OF SP 

The key ingredient in the calculation of stopping power, as outlined in Sec- 
tion II, is the linear response function e{k,uj) of the target material. The 
latter is modelled in most cases by a dense electron gas neutralized by a 
positive background - the so-called jellium model. The effect of lattice is not 
included except perhaps through an effective mass of the electrons. This is a 
reasonable model for target materials like Cu and Al. e{k, cu) for a dense, and 
hence a degenerate electron gas, has been calculated in various approxima- 
tions in the literature. Two of them are: (A) Plasmon-Pole approximation 
(PPA), and (B) the full random phase approximation (RPA). Actually the 
plasmon-pole approximation is a simplification of the RPA response function. 
We shall consider SP in these two approximations. 

3.1 SP in PPA without plasmon dispersion 

Here we consider the simplest model of the dielectric function of a jellium. In 
Ref. [10] (see also the Ref. [15]) a plasmon-pole approximation to e{k,ijj) for 
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an electron gas was used for calculation of dicluster SP. In order to get easily 
obtainable analytical results, Basbas and Ritchie [10] employ a simplified 
form which exhibits collective and single-particle effects 



s{k,uj) 2uj 



H {kc -k)5{uj-ujp) + H{k- kc) S 



2m 



, (10) 



where H{x) is the Heaviside unit-step function, hcUp is the plasma energy of 
the electron gas and the choice = (2mujp/h)^^'^ allows the two 5 functions 
in Eq. (10) to coincide at k = k^ in the k-u plane. 

The first term in Eq. (10) describes the response due to nondispersive 
plasmon excitation in the region k < k^ while the second term describes 
free-electron recoil in the range k > kc (single-particle excitations). This 
approximate function satisfies the sum rule 







Im— -uduj = — ^ (11) 



for all values of k. 

In this approximation if y > {t}jjjp/2mf'^'^ = Vp 

where A ~ V/vp, = l/'n'kpaQ = (4/97r^)"^^^ r^.; = (3/47rnoaQ)"^^^, Uq is 
the electron gas density and = 0.53 x lO^^cm is the Bohr radius, kp is the 
Fermi wave number of the target electrons and Eq = 2.18 GeV/cm. In our 
calculations x {oi r^) serves as a measure of electron density. The result in 
Eq.(12) agrees exactly with the Bethe SP formula, except that the plasmon 
energy of the electron gas hwp appears instead of the usual mean atomic 
excitation energy. Eq. (12) represents the contribution of valence/conduction 
electrons in a solid to the stopping of an ion. 

Using Eq. (10) in Eq. (6), in the high- velocity limit V > Vp {or X'^ > 
x/V^= Aq) one finds 



S,orr{R,^,X) = ^^^i^COs(^-^kpRcOSl9] (13) 
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/■2A dx f \ ( I 

/ — cos\ —kpRcos^d] Jq\ xkpRsm'&x 1 — -— 
Aao ^ V2A ^ ) \ V 4A2 

If one ion trails directly behind the other (t? = 0) from Eq. (13) we find 



Sc,„(fl,O.A) = ^|o„s(3?|..fl)lu(^l+ (14) 

ci {2XkFR) - ci (^-^kpR 

where ci (z) is the integral cosin function 

• / \ f°° 1 cos a; ,^ ^, 

ci(z) = - / dx . (15) 

Jz ^ 

One sees a characteristic oscillatory behavior for large interionic distance 
R. As discussed in [16], fluctuations in the stopping power of a medium for 
a cluster as separation increases are due to electron density variation in the 
wake of the leading ion. The wavelength of these fluctuations is ~ 2t:V /ujp 
for high-velocity projectiles. 

In the case of randomly oriented clusters from Eq. (7) we flnd 



(i?,A) 



where 



Si2 (^^{kpR)^ -Si2{2\{kFR)) 



(16) 



CO 



. , , , sin a; sin(2;) ., , 
si2 {z) = I dx^- = — ^ - ci(^). (17) 

J z ^ 

3.2 SP in PPA with plasmon dispersion 

Plasmons without dispersion are an idealization. In real systems plasmons 
are expected to undergo a dispersion leading to auj(k). The actual dispersion 
( in RPA) can be obtained from the linear response function (see Sec. 3.3). 
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Here we shall utilize a dispersion which is valid for small and intermediate 
values of the wave vector k. Consequently we write 



Im \ = {cu - n{k)) , 



£{k,uj) 2uj 
where the dispersion is given by 



^2/K\ _ , ,2 I ^u2„,2 



Aw? 



In this approximation when 



\b m 

we have, for ISP and CSP, 

So , - 3/5 + V(A2- 3/5)^ -4x73 



(18) 



(19) 



(20) 



'S'ind(A) 



In 



'S'corr(-R, ^1 A) 



7rx2A2 2x/\/3 



r+^^> dx f 
I — 

Jx.(\) X V 2A 



Here kp is the Fermi wave number, and 
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h{x) = 'y/a;4 + —a;2 + 16x73, 



Mx) = W4 (^A2 - x2 - {x^ + 16x73), 



x±(A) 



A2 



4%^ 



(21) 

kpR sin 
(22) 

(23) 
(24) 

(25) 



In the case of randomly oriented clusters we find 

So 



>S'corr(-R, A) 



TTX^A^ 



si2 {kFRx-{X)) — si2 {kFRx-^.{X))] 



(26) 
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3.3 Stopping power in RPA 



Now we will derive the analytical expressions for the SP of a dicluster in a 
fully degenerate {T — 0) electron gas. For this purpose we use the exact 
RPA dielectric response function obtained by Lindhard [4] 

e{z, = 1 + 4 ^) + ^M^^ ^)] ' (27) 

where 



f/+ + l 



v. 



+ 



+ l(C/!-l).n 



f/- +1 



f2{z,u) 



g{l-{u-zf); \u-l\^z^u + l 

0; O^z^u-1, 

0; z u + 1, 

^TTu; ^ z ^ 1 — u 



, (28) 



(29) 



Here, as in Refs. [6, 7, 13, 15], we have introduced the following notations 
z = k/2kF, u — uj/kvpi U± — u ± z. With these notations Eqs. (5), (6) and 
(7) read 



'S'ind(A) 



6So 



oo p\ 
3, 



z-^dz 



f2{z, u)udu 



[z^+xyi{z,u)f + x'fi{z,uy 



(30) 



'5'corr 



6Sn 



/ z^dz 
Jo Jo 



f2{z, u)udu 



vrVA^io "^"io [z^ + x'h{z.u)f + x'fi{z ,u) 

(2^kFRcos^^ Jo (^zkrEsin^^l- ,(31) 



X cos 



^corr(i?, ^) = Sin (2^^i?^) Z^d^ 



/2(^, u)udu 



(32) 

In order to evaluate the integrals by z in Eqs. (30)-(32) at A < 1 (in 
the low-velocity limit) we split the integration region into two domain: ^ 
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z ^ 1 — u and 1 — u^z^l + u, where Im £ ~ /2 7^ 0. However 
at A > 1 (in high-velocity hmit) we need to take into account the region 
l^u^X, O^z^u — 1, where /2 may vanish. The integration in this 
region includes the excitation of collective plasma modes (plasmons) by fast 
charged particles. Consequently, although /2 = the integrals in this region 
are not equal to zero. A calculation of the collective part of SP is facilitated 
if we use the following known expression 



X^f2(z,u) 



[z^ + x'fi{z,u)Y + x'mz,u) 



nd {z' + x'fi{z,u)) = 

5{Z~ Zr{x,u)) 



(33) 



TT 



2 dfi{z,u) 



Z = Zr{x,u) 



where 2;j.(x, m) is the solution of the dispersion equation e{k,uj) = in vari- 
ables z and u. 

Figure 19 shows the solution ^^.(x, «) for various values of x (solid line, 
X — 0.5; dashed hne, x — 0-15; dotted line, x — 0.05). It may be noted 
that the integration domain O^m^A, ^>m + 1, where /2 = 0, does not 
contain the dispersion curve Zr{x,u) calculated for metallic densities x ~ 0.5 
(rs ~ 2). Consequently the SP in this region of variables z and u vanishes 
and there is no plasmon excitation. 

Let us consider the low- velocity limit {V <S vp) of Eqs. (30) and (31). 
In this limit one can obtain simpler expressions for SP. Prom Eqs. (30) and 
(31) we have 



^ind(A)~^A 



z^dz 



^X' Jo [z' + x'm] 



2 ' 



(34) 



S'corr(A, R, i?) 



where 



3Sn , zMz 



:A 



27rx^ Jo [z' + x'f{z)Y 



[^lizkpR) + ^2{zkFR) sin^ §] , 

(35) 

(36) 



e'-- sin (20+ e cos (20 
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(37) 



(38) 



From Eqs. (36) and (37) it follows that <l>i(0 ^ 2/3, $2(0 at ^ 
and consequently, as expected, 5'corr(A, i?, •«?) — >■ 'S'ind(A) when R ^ 0. Note 
that (as is well-known [4, 13]) in the low- velocity limit the SP is proportional 
to the velocity of particle (Eqs. (34) and (35)). Thus the vicinage function 
r(A, i?, t?) at A <^ 1 depends only on interionic distance R and orientation 
angle 'd. 



4 NUMERICAL RESULTS AND DISCUS- 



Using the theoretical results of Sees. 2 and 3, we have made extensive nu- 
merical calculations of stopping power (SP) and related quantities. In this 
section we present detailed numerical results for two target materials, Al 
and Cu. These two targets have been chosen because of their frequent use 
in experiments and also because of their different electron densities. In our 
calculations x (or ^g) is a measure of electron density. 

As a simple but generic example of a projectile, we have considered 
a diproton cluster for which we present theoretical results for the follow- 
ing quantities of physical interest: stopping power (SP/2), vicinage func- 
tion (VP), angle-averaged stopping power (ASP/2), angle-averaged vicinage 
function (AVF) together with the dependence of SP/2 and VP on R, the 
inter-ionic separation distance within the cluster. The reason why SP has 
been divided by a factor of 2 is that the SP results for a diproton cluster 
are expected to reduce asymptotically (as R tends to infinity) to those for 
two uncorrelated protons, the latter being referred to as ISP. ASP has been 
treated in the same way. 

In our calculations of these quantities we have employed the linear re- 
sponse approach which assumes a swift ion-cluster projectile and also that 
the ion cluster presents a weak perturbation on the target plasma. The va- 
lidity of the linear response approach to study ion-cluster stopping has been 
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discussed in detail by Zwicknagel and Deutsch [17]. We refer the reader to 
their insightful discussion. 

We model the Al and Cu targets by a dense (degenerate) electron gas neu- 
tralized by a positive background (the jellium model) with electron densities 
appropriate for the respective targets. The hnear response of the target elec- 
tron gas, which couples the cluster projectile to the target, is considered at 
three levels of approximations to the dielectric function e{k, u) as discussed in 
sections 2 and 3. In the context of stopping power these approximations are 
subject to the following general remarks: The plasmon-pole approximation 
(PPA) is valid only in the high velocity regime when the mean velocity V of 
the cluster is > vp, the Fermi velocity of the target electrons. For V < vp 
and for velocities near the threshold of collective mode excitations, this ap- 
proximation is not adequate. The RPA overcomes this limitation although it 
cannot account for short-range correlations in the electron gas. Within PPA 
itself, PPA-1 (without plasmon dispersion) is more limited than PPA-2 (with 
plasmon dispersion). The figures we present serve as a comparative study 
of how these levels of approximation affect the various physical quantities 
related to stopping power. 

Figs. 1-4 show cluster stopping power (CSP and its dependence on var- 
ious quantities of experimental interest) . Let us first note that these figures 
are presented for two specific values, and 7r/2, of the angle Correlations 
between the two ions in the dicluster are maxiniTim and minimum, respec- 
tively, for these two values of i}. The objective is then to see how, for these 
maximum and minimum configurations, CSP depends on R and V/vp. Fig. 
1 shows CSP for Al target with R — 10~^cm, as a function oiV/vp for the 
two above-mentioned values of i?, within PPA. The lines without circles cor- 
respond to PPA-1 and, with circles, to PPA-2. The angular dependence of 
CSP is particularly noteworthy. It is seen that in a medium velocity range 
{V < 2vf), CSP has a remarkably higher value for the larger value oii). This 
is likely due to single-particle excitations in this velocity range. In the higher 
velocity range, the dicluster wake-field excitations become important and we 
find that the situation is reverse in the higher velocity range {V > 2vf) for 
which CSP for -i? = is larger than for '& = 7t/2. 

In the low velocity range the difference between PPA-1 and PPA-2 (for 
both d — and 7r/2) is noticeable while in the high velocity range this differ- 
ence becomes negligible. This is again due to single-particle excitations in the 
low velocity range. For comparison, we have also presented the uncorrelated 
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stopping power (ISP). 

When we increase the inter-ionic separation distance R from 10^*^cm to 
5 X 10~^cm, keeping other physical parameters the same, some interesting 
changes occur, as can be seen from Fig. 2. A noticeable change is that 
now, for y < 2vf, CSP for i? = is higher than that for i? = 7r/2. This 
sensitivity of CSP to the angle -& as Ris varied may be due a combination of 
factors. The dicluster behaves like a compact project for small R, and like 
an extended projectile for large R. This has a bearing on S'corr given in Eqs. 
(13) and (22). Correlation effects are expected to be maximum when the 
two ions are aligned with each other in the direction of propagation of the 
dicluster projectile motion (i? = 0) while they decay (at least for V > vp) 
when i) tends to n/2, the latter behavior being related to the wake-field due 
to the leading ion. The oscillation amplitude in Scon- tends to decrease from 
i) = to = 7r/2 (the Cherenkov cone). However when R is small each 
ion is influenced by the unscreened held of the other ion. For model solid 
targets, the Cherenkov cone semivertex is = arcsin(-\/(l6i>^/y) [18]. '&c 
approximately equal to 22.8°, 7.4°, and 0.08° for V = 2vf, V = Gvp, and 
V = lOvp, respectively. Consequently in the high velocity range the trailing 
ion moves inside the Cherenkov cone of the leading ion only for almost aligned 
diclusters. The behavior of CSP shown in Figs. 1 and 2 reflects these features 
within the hnear response and for PPA. It will be noted that the high values 
of SP are due to the PPA-1 approximation. PPA-2 decreases these values 
to a small extent. Later, when we use a more realistic, namely RPA, for 
the linear response function (Figs. 20 and 21) SP considerably decreases in 
strength. 

Figs. 3 and 4 show SP for Cu, another commonly used metallic target. 
These figures show patterns similar to those in Figs. 1 and 2, except that 
CSP and ISP have lower values over the entire range oiV/vp- This is because 
Al has a higher electron density than Cu. 

The vicinage function (VF) given by Eq. (8) has been plotted as a func- 
tion of the beam velocity for Al target in Figs. 5 and 6. This function shows 
an interplay between i) and R more strikingly than CSP. Figs. 7 and 8 display 
a similar behavior for Cu target. 

As stated in Sec. 2, an average stopping power (ASP) is of experimental 
interest. Figs. 9 and 10 show ASP for Al and Cu, respectively. ISP is also 
shown, for comparison. The role of PPA-1 and PPA-2 is now more clearly 
seen. 
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In the same spirit we have plotted AVF for Al and Cu in Figs. 11 and 

12. The role of R is highlighted in these figures. However it will be noticed 
that PPA-1 and PPA-2 make practically no distinction for AVF. 

We have so far plotted SP or ASP (divided by a factor of 2 in both the 
cases) vs the beam velocity V/vp, for some values of the separation distance 
R. We now look for some complementary information about SP, and plot SP 
as a function of R with V = 3?,'^,% for Al target. Fig. 13 shows an oscillatory 
character of SP with respect to R. The oscillations are the highest for i) = 
and lowest for = 7r/2. The role of PPA-1 and PPA-2 is clearly seen for 
-& — 0. Fig. 14 shows a similar behavior of SP for Cu although the amplitudes 
are now weaker. 

In the same way the vicinage function (VF) is plotted in Figs. 15 and 16, 
for Al and Cu targets, respectively. 

Fig. 17 shows ASP vs R for both Al and Cu targets. The difference 
between PPA-1 and PPA-2 is negligible and the Cu target has ASP smaller 
than for Al. Now, there is something interesting about Fig. 18 which shows 
AVF. The difference between PPA-1 and PPA-2 is again negligible. But 
let us note that data for both Al and Cu lie practically on the same curve! 
Recalling the definition of VF, Eq. (8) one can see from Eqs. (12), (16), (21) 
and (26) that AVF has a weak dependence on target density. Also, when 
A = V/vf > 2, jS'ind does not noticeably depend on PPA-1 and PPA-2. These 
features combine to lead to the behavior of AVF as seen in Fig. 18. 

We have so far presented results for PPA-1 and PPA-2. A more realis- 
tic linear response function, namely the exact random phase approximation 
(RPA) will now be used for the metallic target. The theoretical results for 
SP etc. have been presented in Sec. 3.3. As part of calculating SP in RPA 
it is useful to examine the plasmon dispersion obtained through e{zr, u) = 0, 
where z and u have been defined in Sec. 3.3. Fig. 19 displays Zr{x,u) vs u 
for three electron density parameter values. 

Next we present ISP and CSP in RPA, for Al and Cu in Figs. 20 and 
21, corresponding to R — 10~^cm and R — 5x lO^^cm, respectively. For the 
sake of a better presentation of the data we have also separately displayed 
the Al data in Figs. 20a and 21a. The RPA results show that SP and ISP 
decrease in strength with an improved linear response function. This should 
be of relevance to experiments. 

Next, VF in RPA vs V/vp is presented in Fig. 22, for Al (fines without 
circles) and for Cu (fines with circles), corresponding to R — 10~^cm and 
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for 1? = and 7r/2. This figure may be compared with Figs. 5 and 7. For 

V/vp < 2, the curves for VF in RPA tend toward finite values whereas the 
VF-curves in PPA do not although the angular trend is similar. A similar 
contrast may be noted between Figs. 6 and 8, and Fig. 23, corresponding to 
R — bx 10~^cm. Again, these findings are of experimental relevance. 

Averaged SP vs. V/vp in RPA is presented in Fig. 24, for Al (curves 
without circles) and for Cu (curves with circles) along with ISP, correspond- 
ing to R = 10~^cm and 5 x 10~®cm. This figure may be compared with Figs. 
9 and 10. Fig. 24 shows an expected overall decrease in the strength of ASP 
in RPA. 

There are similarities but also some interesting differences if we compare 
Figs. 11 and 12 with Fig. 25. The latter shows AVF in RPA for Al (curves 
without circles) and for Cu (curves with circles). The differences are more 
noteworthy for R = 10~*cm. 

SP vs R in RPA is plotted in Fig. 26, corresponding to V — Svp, for Al 
and Cu in the previously stated scheme. When Fig. 26 is compared with 
Figs. 13 and 14, the differences between PPA and RPA become particularly 
striking. 

A similar contrast is provided by a comparison of Fig. 27 with Figs. 15 
and 16, for VF vs. R in RPA and PPA, corresponding to V — Svp and for 
■J? = and 7r/2. For comparison AVF is also plotted in Fig. 27. 

This completes our extensive presentation of figures exhibiting various 
aspects of the stopping power of a diproton cluster in PPA and RPA, for Al 
and Cu targets. 

5 SUMMARY 

In this paper we have presented a comprehensive theoretical study of stopping 
power (SP) of a dicluster of protons in a metallic target. After a general 
introduction to SP of a cluster of two point-like ions, in Sec. 2, theoretical 
calculations of SP based on the linear response theory and using PPA without 
and with plasmon dispersion and then with RPA are discussed in Sec. 3. The 
theoretical expressions for a number of physical quantities derived in section 
lead to a detailed presentation, in Sec. 4, of a large collection of data through 
figures on correlated stopping power (CSP), vicinage function (VF), average 
stopping power (ASP) and average vicinage function (AVF) of a diproton 
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cluster projectile for two metallic targets, Al and Cu. Whenever relevant, 
wc have also provided a plot of independent (i.e. single-ion) stopping power 
(ISP) for comparison. 

With the proviso stated in Sec. 4, SP and related quantities have been 
studied within a linear response formalism; some analytical and all numerical 
results have been obtained corresponding to three approximations to the di- 
electric function of the target electron gas-the plasmon-pole approximation 
(PPA) without dispersion (PPA-1) and with dispersion (PPA-2), and also 
with the random phase approximation (RPA). To our knowledge this is the 
most comprehensive calculation of the SP-related physical quantities using 
all the three dominant approximations to the linear response function. The 
results we have presented demonstrate that with regard to several physical 
quantities of primary interest the difference between PPA and RPA is sub- 
stantial while for others, specially for average quantities, this difference may 
not be of practical significance. 

It will be of interest to go beyond RPA in order to include some short- 
range correlations in the electron gas and to study how dicluster SP is af- 
fected. However calculating the linear response function by including electron 
energy bands is rather involved and detailed theoretical studies of SP with 
band structure effects included have not yet been reported in the literature. 
One can include some aspect of band structure in a rather approximate man- 
ner through an effective mass for the electrons. 

Another aspect we have not considered in this paper is some effect of 
disorder in the target medium. In real metals electrons suffer collisions with 
impurities etc. We intend to address this issue in the context of stopping 
power in a separate study. 
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Figure Captions 

Fig. 1. SP/2 of a diproton cluster with R = 10~^cm, vs V/vp for Al 
target (r^ = 2.07). = (dotted line), 1^ = 7i/2 (dashed line); ISP (solid 
line) . The lines with and without circles correspond to PPA with and without 
dispersion, respectively. 

Fig. 2. SP/2 of a diproton cluster with R = 5 x 10~®cm, vs V/vp for Al 
target, i) = (dotted line), i) = 7i/2 (dashed line); ISP (solid line). The lines 
with and without circles correspond to PPA with and without dispersion, 
respectively. 

Fig. 3. SP/2 of a diproton cluster with R = 10~^cm, vs V/vp for Cu 
target (r^ = 2.68). 19 = (dotted hue), 13 = 7r/2 (dashed line); ISP (solid 
line). The lines with and without circles correspond to PPA with and without 
dispersion, respectively. 

Fig. 4. SP/2 of a diproton cluster with i? = 5 x 10~^cm, vs V/vp for Cu 
target. -& — (dotted hne), ■»? = 7r/2 (dashed line); ISP (solid hne). The hues 
with and without circles correspond to PPA with and without dispersion, 
respectively. 

Fig. 5. VF of a diproton cluster with R = 10~^cm, vs V/vp for Al target. 
-& — (solid line), = 7r/2 (dotted line). The hues with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 6. VF of a diproton cluster with R = 5 x 10~^cm, vs V/vp for Al 
target. i9 = (solid line), d = ti /2 (dotted line). The lines with and without 
circles correspond to PPA with and without dispersion, respectively. 

Fig. 7. VF of a diproton cluster with R — 10~^cm, vs V/vp for Cu target. 
■J? = (solid line), •& — (dotted line). The hues with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 8. VF of a diproton cluster with R = 5 x 10~^cm, vs V/vp for Cu 
target. i3 = (solid line), i9 = 7t/2 (dotted line). The lines with and without 
circles correspond to PPA with and without dispersion, respectively. 

Fig. 9. ASP/2 of a diproton cluster with R = 10~^cm (dotted hne) 
and R = 5 X 10~^cm (dashed line) vs V/vp for Al target. Solid line, ISP. 
The lines with and without circles correspond to PPA with and without 
dispersion, respectively. 

Fig. 10. ASP/2 of a diproton cluster with R = 10-^cm (dotted hne) 
and R — 5 X 10~^cm (dashed hne) vs V/vp for Cu target. Solid line, ISP. 
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The lines with and without circles correspond to PPA with and without 
dispersion, respectively. 

Fig. 11. AVF of a diproton cluster with R = lO^^cm (solid line) and 
R — 5 X 10~^cm (dotted line) vs V/vp for Al target. The lines with and 
without circles correspond to PPA with and without dispersion, respectively. 

Fig. 12. AVF of a diproton cluster with R — 10~^cm (solid hue) and 
R = 5 X 10~^cm (dotted line) vs V/v^ for Cu target. The lines with and 
without circles correspond to PPA with and without dispersion, respectively. 

Fig. 13. SP/2 of a diproton cluster with V = Svp vs R for Al target. 
-& — (solid line), -d = n/2 (dotted line). The hues with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 14. SP/2 of a diproton cluster with V = 3vf vs R for Cu target. 
= (solid line), d = ti/2 (dotted line). The lines with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 15. VF of a diproton cluster with V — ?>vf vs it! for Al target. 
■J? = (solid line), ■& — 'njl (dotted line). The hues with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 16. VF of a diproton cluster with V = Svp vs R for Cu target. 
"d = (solid line), = 7r/2 (dotted line). The lines with and without circles 
correspond to PPA with and without dispersion, respectively. 

Fig. 17. ASP/2 of a diproton cluster with V = 3^^ vs R for Al (the hues 
with square symbols) and Cu (the lines with circles) targets. Dotted and 
solid lines, PPA with and without dispersion, respectively. 

Fig. 18. AVF of a diproton cluster with V = 3vp vs R for Al (the lines 
with square symbols) and Cu (the lines with circles) targets. Dotted and 
solid lines, PPA with and without dispersion, respectively. 

Fig. 19. Relation between Zr{X: and u, as obtained from the dispersion 
equation e {zr, u) — 0. Solid line: x — 0-5, dashed line: x — 0.15, dotted line: 
X = 0.05. 

Fig. 20. SP/2 of a diproton cluster with R = lO'^cm in a RPA vs V/vp 
for Al (the hues without symbols) and Cu (the lines with circles) targets. 
d = (dotted line), i) = n/2 (dashed line); ISP (solid line). 

Fig. 20a. SP/2 of a diproton cluster with R = lO^^cm in a RPA vs V/vp 
for Al target, = (dotted line), -d = ti/2 (dashed line); ISP (solid line). 

Fig. 21. SP/2 of a diproton cluster with R = 5 x 10~^cm in a RPA 
vs V/vp for Al (the hues without symbols) and Cu (the hues with circles) 
targets. — (dotted line), -& — 7r/2 (dashed line); ISP (sohd line). 
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Fig. 21a. SP/2 of a diproton cluster with i? = 5 x lO'^cm in a RPA vs 
V/vf for Al target. -(9 = (dotted line), ^ = 7r/2 (dashed line); ISP (solid 
line) . 

Fig. 22. VF of a diproton cluster with R = 10~®cm in a RPA vs V/vp for 
Al (the hues without symbols) and Cu (the lines with circles) targets. — 
(solid hne), -& — 7r/2 (dotted hne). 

Fig. 23. VF of a diproton cluster with i? = 5 x 10~*cm in a RPA vs V/vp 
for Al (the lines without symbols) and Cu (the lines with circles) targets. 
= (solid line), = 7r/2 (dotted line). 

Fig. 24. ASP/2 of a diproton cluster with R = 10~^cm (dotted hne) and 
R — 5 X 10~^cm (dashed line) in a RPA vs V/vp for Al (the hues without 
symbols) and Cu (the lines with circles) targets; ISP (solid line). 

Fig. 25. AVF of a diproton cluster with R = 10"^cm (solid line) and 
R — 5 X 10~^cm (dotted line) in a RPA vs V/vp for Al (the lines without 
symbols) and Cu (the lines with circles) targets. 

Fig. 26. SP/2 of a diproton cluster with V — 3vp in a, RPA vs R for 
Al (the lines without symbols) and Cu (the lines with circles) targets; ASP 
(solid line), = (dotted hne), d = ti /2 (dashed line). 

Fig. 27. VF of a diproton cluster with V = "ivp in a RPA vs R for 
Al (the hues without symbols) and Cu (the hues with circles) targets; AVF 
(solid hne), ■& — Q (dotted hne), ■& — (dashed line). 
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